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Abstract 

We establish new, optimal gradient continuity estimates for solutions to a class of 2nd order 
partial differential equations, L{X,Vii,D^u) — /, whose diffusion properties (ellipticity) degen- 
erate along the a priori unknown singular set of an existing solution, S(m) := {X : Vu{X) = 0}. 
The innovative feature of our main result concerns its optimality - the sharp, encoded smooth- 
ness effects of the operator Such a quantitative information usually plays a decisive role in the 
analysis of a number of analytic and geometric problems. 
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1 Introduction 

Regularity theory for solutions to partial differential equations has been a central subject of research 
since the foundation of the modem analysis of PDEs, back in the 18th century. Of particular interest 
are physical and social phenomena that involve diffusion processes, whose mathematical models are 
governed by second order elliptic PDEs. 

Smoothness of weak solutions to 2nd order uniformly elliptic equations, both in divergence and 
in non-divergence forms, is nowadays fairly well established. The cornerstone of the theory is a 
universal modulus of continuity for solutions to the homogeneous equation: Lu = 0. This is the 
contents of DeGiorgi-Nash-Moser theory for the divergence equations and Krylov-Safonov Hamack 
inequality for non-divergence operators. 

Despite of the structurally profound importance of the supra-cited works, a large number of math- 
ematical models involve operators whose ellipticity degenerates along an a priori unknown region, 
that might depend on the solution itself: the free boundary of the problem. This fact impels less 
efficient diffusion features for the model near such a region and therefore the regularity theory for 
solutions to such equations become more sophisticated from the mathematical view point. 

The most typical case of elliptic degeneracy occurs along the singular- set of an existing solution: 

S(m) ■={X ■.Vu{X)=0}. 
In fact, a number of degenerate elliptic equations has its degree of degeneracy comparable to 

(1.1) f{Vu)\D\\^l, 

for some function / : M'' — > M, with Zero (/) = {0}. Thus, understanding the precise effect on the lack 
of smoothness impelled by the model equation dl.ll ) shades lights on the underlying sharp regularity 
theory for a number of typical degenerate elliptic operators. 

The main goal of this present work is to derive sharp, interior regularity estimates for degenerate 
elliptic equations of the general form 

(1.2) J^{X,Vu)F{X,D'^u)=f{X), BiCR'^, 

where / € L°°{B\) and J^: Bi xR'^ degenerates within a comparable rate to a power of the 
magnitude of the gradient, i.e., 

(1.3) A|p|^<^(X,p) <A|p|^ 

for some 7 > 0. The 2"'' order operator F: Bi x Sym(<i) ^ M in equation (11.21 ) is responsible for 
diffusion, i.e., F will be assumed to be a generic fully nonlinear unifoiTnly elliptic operator: Xld^xd < 
dijF{X,M) < AXldcixd- In Section|2]we give a more appropriate notion of ellipticity. 

Regularity theory for viscosity solutions to fully nonlinear uniformly elliptic equations. 



f (D^m) =0, 
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has attracted the attention of the mathematical community for the last thi^ee decades or so. It is 
well established that solutions to the homogeneous equation is locally of class C^'"° for a universal 
exponent «()» depending only on d, X and A, see for instance ||2]. If no additional structure is 
imposed on F, C^ "^ is in fact optimal, see HH, lUl, ll20l . 

A quick inference on the stmcture of equation (ll.2l i reveals that no universal regularity theory 
for such equation could go beyond C^'"". In fact the degeneracy term J^(X, Vi<) forces solutions to 
be less regular than solutions to the uniformly elliptic problem near its singular set. This particular 
feature indicates that obtaining sharp regularity estimates for solutions to (11.21 ) should not follow from 
perturbation techniques. Indeed, it requires new ideas involving an interplay balance between the 
universal regularity theory for uniform elliptic equations and the degeneracy effect on the diffusion 
attributes of the operator compelled by (ll.Bl l. 

In this present work we show that a viscosity solution, u, to (11.21 ) is pointwise differentiable and its 
gradient, Vm, is locally of class (7°'™"{"o 1+7} jhe sharp Holder continuity exponent for the gradient 
of solutions, 

(1.4) i3:=min{ao-,y^}, 

gives the precise optimal and universal regularity for degenerate equations of the type (ll.21 i. We 
further estimate the corresponding maximum regularity norm of w by a constant that depends only 
on universal parameters, 7, ||/||oo and ||m||oo. We have postponed the precise statement of the main 
Theorem to Section |3] 

The key, innovative feature of our main result lies precisely in the optimality of the gradient 
Holder continuity exponent of a solution to the degenerate equation (11.21 ). which in turn is an impor- 
tant piece of information in a number of qualitative analysis of PDEs, such as blow-up analysis, free 
boundary problems, geometric estimates, etc. It is quantitative bonus acquisition to the recent result 
of Imbert and Silvestre, lITTl . where it is proven that viscosity solutions to (11.21 ) are continuously 
differentiable. The logistic reasoning of the proof of our main result is inspired by recent works of 
the third author, Il23ll . Il24l . Il25l . and it further uses the main crack from lITTI to access a priori 
estimate for solutions to (ll.21 i. 

The paper is organized as follows. In Section |2] we gather the most relevant notations and known 
results we shall use in the paper. In Section |3] we present the main Theorem proven in this work. 
In Section |4] we provide a few implications the sharp estimates from Section |3] have towards the 
solvability of some well known open problems in the elliptic regularity theory. The proof of Theorem 
13. H is delivered in the remaining Sections |5]|6] Eland [8] 

2 Notation and preliminares 

In this article we use standard notation from classical literature. The equations and problems studied 
in this paper are modeld in the J-dimensional Euclidean space, . The open ball of radius r > 
centered at the point Xq is denoted by BriXo). Usually ball of radius r, centered at the origin is written 
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simply as B,-. For a function m: Si — )• R, we denote its gradient and its Hessian at a point X G Bi 
respectively by 

Vu{X) := {dju)i<j<d and D^u{X) := {diju)i<ij<d, 

where dju and djjU denote the j'-th directional derivative of u and the i-th directional derivative of 
dju, respectively. 

The space of all dxd symmetric matrices is denoted by Sym(<i). An operator F : BiX Sym{d) — )■ 
M is said to be uniformly elliptic if there there exist two positive constants < A < A such that, for 
any M G Sym{d) and X eBi, 

(2.1) X\\P\\ <F{X,M + P)-F{X,M) <A\\P\\, VP > 0. 

Any operator F satisfying the ellipticity condition (12.11 ) will be referred hereafter in this paper as 
a (A, A) -elliptic operator. Also, following classical terminology, any constant or entity that depends 
only on dimension and the ellipticity parameters A and A will be called universal. For normalization 
purposes, we assume, with no loss of generality, throughout the text that F(X,0) = 0, VX G B^. 

For an operator G: Bi xW^ x Sym(J) R, we say a function u G C*'(Bi) is a viscosity super- 
solution to G{X,Vu,D'^u) = 0, if whenever we touch the graph of u by below at a point F G Si by 
a smooth function <p, there holds G{Y ,'V (p{Y) , (p{Y)) < 0. We say u G C"(Bi) is a viscosity sub- 
solution to G{X,Vu,D^u) = 0, if whenever we touch the graph of u by above at a point Z G Bi by a 
smooth function 0, there holds G{Y,'V(p{Y),D^(p{Y)) > 0. We say w is a viscosity solution if it is a 
viscosity super-solution and a viscosity sub-solution. The crucial observation on the above definition 
is that if G is non-decreasing on M with respect to the partial order of symmetric matrices, then 
the classical notion of solution, sub-solution and super-solution is equivalent to the corresponding 
viscosity terms, provided the function is of class C^. The theory of viscosity solutions to non-linear 
2nd order PDEs is nowadays fairly well established. We refer the readers to the classical article ||6]. 

Let us discuss now a little bit further the existing regularity theory for uniformly elliptic equations. 
As mentioned earlier in the Introduction, it follows from the celebrated Krylov-Safonov Harnack 
inequality, see for instance lUl, that any viscosity solution to the constant coefficient, homogeneous 
equation 

(2.2) F{D^h)=0, 

is locally of class C^'"" for a universal exponent < tto < 1- Hereafter in this paper, Oq = ao{d , X , A) 
will always denote the optimal Holder continuity exponent for solutions constant coefficients, ho- 
mogeneous, (A, A) -elliptic equation (12.21 ). If no extra structural condition is imposed on F, C^^^' 
is indeed the optimal regularity possible, lITSl . |[T9l . EOl . However, under convexity or concavity 
assumption on F, solutions are of class C^. This is a celebrated result due to Evans |7J and Krylov 
|[15l[T6l, independendy. 

For varying coefficient equations, solutions are in general only c{)^^' and this is the optimal 
regularity available, unless we impose some continuity assumption on the coefficients, i.e., on the 
map Z I— )• f (X, •), Such condition is quite natural and it is present even in the linear theory: Lu := 
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aij{X)DijU. Since we aim for a universal C^ " estimate for solutions to equation (11.21 ). hereafter we 
shall assume a uniform continuity assumption on the coefficients of F, appearing in (11.21 ). namely 

\F(X,M)-F(Y,M)\ 
(2-3) sup ' ^ ^ ^' <C(0i\X-Y\), 

\\M\\<1 II™ II 

where C > is a positive constant and O) is a normalized modulus of continuity, i.e., CO : 1R+ — M+ is 
increasing, (o(0+) = and co(l) = 1. Such condition could be relaxed: it suffices some sort of VMO 
condition, see |[T]|. We have decided to present the results of this present article under (|2.3l l for sake 
of simplicity. For notation convenience, we will call 

(2.4) := inf |c> : sup '^^^'^h <C(oi\X -Y\), yXJ e bA . 

[ ||M||<1 ll^ll J 

We close this Section by mentioning that under continuity condition on the coefficients, viscosity 
solutions to 

F{X,D^u)=f{X)eL-{Bi) 

are locally of class C^'^ , for any < j3 < OCq, where (Xo is the optimal Holder exponent for solutions 
to constant coefficient, homogeneous equation F{D^h) = 0, coming from Kiylov-Safonov, Caffarelli 
universal regularity theory. See |jl], [|24l . 



3 Main results 

In this Section, we shall present the main result we will prove in this present work. As mentioned 
earUer, the principal, ultimate goal of this article is to understand the sharp smoothness estimates for 
functions u, satisfying 

(3.1) |VM|^-|F(D2ii)|<l, 0>O, 

in viscosity sense, for some uniformly elliptic operator F . Clearly, as commented in the previous 
Section, even in the non-degenerate case, 0=0, the best regularity possible is C^^^ " . The delicate 
point, though, is to obtain a universal estimate, fine enough as to sense and deem the singularity 
appearing in RHS of Equation ( 13.11 ). along the singular set (Vm)^^(O), as > varies. 
As to grasp some feelings on what one should expect, let us naively look at the ODE 

u"{t) = {u')-^, m(0) =m'(0) =0, 

2+e 

which can be simply solved for t G (0, oo). The solution is u{t) = ? i+e . After some heuristics inference, 

2+9 

it becomes reasonable to accept that C i+e is another upper barrier for any universal regularity estimate 
for Equation (13.11 ). Thus, if no further obscure complexity interferes on the elliptic regularity theory 
for fully nonlinear- degenerate elliptic equation, the ideal optimal regularity estimate one should hope 
for functions satisfying Equation ( 13.11 ) should be c''™"^"o T+e}. 

After these technical free and didactical comments, we are able to state the main result we estab- 
lish in this paper. 
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Theorem 3.1. Let ube a viscosity solutions to 



(3.2) 



J^{X,Vu)F{X,D^u)=f{X) in Bi. 



Assume f ^U°{B\), satisfies (11.31 ) and F : Bi x Sym{d) — )• M is uniform elliptic with continuous 
coefficients, i.e., satisfying (12.11 ) and (12.31 ). Fixed an exponent 



there exists a constant C((i,A,A,7, ||f ||/||c,o,o;) > 0, depending only on t/, A,A, 7, ||f Hm, ||/|| 
and a, such that 

||m|Ici-«{b,/2) < C((i,A,A,7, ||f lloj, ||/||oo,a) • 
An important consequence of Theorem 13. H is the following: 
Corollary 3.2. Let u be a viscosity solutions to 



Assume f € L°°{B\ ), M' satisfies (11.31) . F is uniformly elliptic and concave. Then u is locally in C ' 
and this regularity is optimal. 

Corollary l3.2l follows from Theorem 13. 1 1 since solutions to concave equations are locally of class 
C''^ by Evans-Krylov Theorem. 

It is interesting to understand Theorem 13. H as a model classification for degenerate elliptic equa- 
tions, linking the magnitude of the degeneracy of the operator to the optimal regularity of solutions. 
A quantitative, intrinsic signature of the degeneracy properties of the equation. In fact, as mentioned 
earlier, many classical equations have their degree of degeneracy comparable to a model equation of 
the form |VM|^|f (D^m)| < 1. We shall explore this perspective within the next Section. 

4 Applications and further insights 

The heuristics from the "degeneracy classification" mentioned in the previous paragraph has indeed 
a wide range of applicability. In this intermediary Section we comment on some consequences the 
optimal regularity estimates stated in Section |3] have in the elliptic regularity theory. 

In the sequel we shall use Theorem 13. II and Corollarv l3.2| to solve particular cases of some well 
known open problems. The results provided in this section give hope that decisive progress can be 
attempted for the general cases in the near future. 

4.1 Equations from the theory of superconductivity 

We start off by commenting on some applications Theorem l3.1l has to the theory of superconductivity, 
where fully nonlinear equations with patches of zero gradient 




(3.3) 



J^{X,Vu)F{D^u) = f{X) in Bi. 



(4.1) 



F{X,D\i)=g{X,u)X{\vr,\>o] 
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governs the mathematical models. Equation (14.11 ) represents the stationary equation for the mean field 
theory of superconducting vortices when the scalar stream function admits a functional dependence 
on the scalar magnetic potential, see fS]. Existence and regularity properties of Equation (14. H were 
studied in [3] and in 141. The novelty to study Equation (14.1b is that one tests the equation only 
for touching polynomials for which |VP(Xo)| 7^ 0. It is proven in Q, Corollary 7, that solutions 
are locally C°'" for some < a < 1. For concave operators, it is proven, see ||3] Corollary 8, that 
solutions are in W^ '^. An application of Alt-Caffarelli-Friedman monotonicity formula, ||3] Lemma 
9, gives regularity C' ' for the particular problem 

(4.2) Au = cu X{\vu\>0}- 

Equation (14.11) can be obtained as a limiting problem, as 5 ^ 0, for the family of singular equa- 
tions 

(Es) \Vusf-FiX,D\s) = g{X,us), Si. 

Indeed, it follows from Theorem [3.1 1 that if u§ is a normalized solution to ( [Eg] ), for d small enough, 
i.e., for 



5 < 1 - , 



then we can estimate 



(4-3) ||m|| Ua, <C, 



for a constant C > 1, that does not depend on 5. In particular, estimate (14.31 ) gives local compactness 
for the family of solutions {ug}syQ to ^Eg^. Let mq be a limiting point of such a sequence, i.e. 



Mo = limMg , 



for 5y = 0(1). From (14.31 ). we have. 



(4.4) VMg^ — > Vmq locally uniformly, 

(4.5) MO € c/„f (Si). 
Now, fixed a regular point Z G Si of uq, i.e., 

|VMo(Z)j >0. 

Gradient convergence (14.41 ) and estimate (14.51 ) yield the existence of a small T] > 0, such that 



inf \Vus \ > -^|Vmo(Z)| =: cq, 
B^(Z) 10 

for all 5 < 1. Thus, 

g{X,us) -iVusr^ — >g{X,uo), uniformly in S^(Z). 
We sum up the above discussion as the following Theorem: 
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Theorem 4.1. Let ug G C^{Bi) be a viscosity solution to ^Eg^ , with < 1, 5 ^ 1, where g is 
continuous w.r.t. u and measurable bounded w.r.t. X. Assume the operator F is under the hypotheses 
ofTheorem \3.1\ Then, fixed a number a < (XQ{d,X,A), for 5 small enough, we have 

ll"5llci-«(B4/5) <C{d,X,A,a). 

In particular, 

and uq is a viscosity solution to (14. It . 

The advantage of Theorem 14.11 in comparison to the regularity theory developed in ||3l is that it 
provides the asymptotically sharp C' " estimate in the general case of fully nonlinear operators, not 
necessarily concave. 

4.2 Visiting the theory of oo-laplacian 

Let us now visit the theory of the oo-laplacian operator, i.e., 

(4.6) AooV := 52v,VyV;y 

'J 

which is related to the problem of best Lipschitz extension to a given boundary datum - a nonlinear 
and highly degenerate elliptic operator. The theory of infinity-harmonic functions, i.e., solutions to 
the homogeneous PDE 

AJi = 0, 

has received a great deal of attention. One of the main open problems in the modem theory of PDEs is 
whether infinity-harmonic functions are of class C' . This conjecture has been answered positively by 
O. Savin |[22l in the plane. Evans and Savin, HI sharpened the result to C^ " for some small a > 0, but 
still in only dimension two. Quite recently, Evans and Smart proved that infinity-harmonic functions 
are everywhere differentiable regardless the dimension, iQ. Nevertheless, no continuity feature of Vu 
can be inferred by their ingenious arguments. The famous example of the infinity-harmonic function 

(4.7) a{x,y):=x-3 -y-^ 

due to Aronsson from the late 60's sets the ideal optimal regularity theory for such a problem. That 
is, no universal regularity theory for infinity harmonic functions can go beyond C''3. Up to our 
knowledge, there has been no prior meaningful mathematical indication that infinity-harmonic func- 
tions should or should not have a universal C^ a regularity theory, other than speculation based on 
Aronsson's example. Another way, though, to surmise the C''5 conjecture for infinity harmonic 
function would be by exploring the scaling properties of the equation. For instance, it one writes the 
infinity-laplacian as 

A^v = {Vvy -D^vVv, 
it becomes tempting to compare its degeneracy feature with 

(4.8) |VMp.|AM|<l, 
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1 _±_ 

that has the same scaling properties as Aoo and whose solutions are locally C ' 1+2 regular from Corol- 



lary |32] Although, it is not in general true that infinity-harmonic functions satisfy (I4.8I ). this obser- 
vation sets an interesting heuristic guide. 

Notice that Aronsson's example - as many popular examples in the theory of PDEs - is a function 
of separable variables. In the sequel, as an application of Corollary 13.21 we show that any infinity- 
hamionic function with separable variables is locally of class 

Proposition 4.2. Let ir. B\ C M'' — ?• M Z^e infinity harmonic. Assume u is a fiinction of separable 
variables, i.e., 

U{X) = Oi{xi) + 02{x2) + ---Od{Xd), 

for Oi G C"(Bi). Then u G C''3 {3^/2). 
Proof. Formal direct computation gives 

(4.9) = AooM = |ai(xi)f af(^i) + |ai(x2)f a^'(x2) + • • • + \ad{xd)'\^o'Ji{xd). 

It is a manner of routine to justify the above computation using the viscosity solution machinery. We 
notice, however, that the /th term in (|4.91 l depends only upon the variable Xj. Thus, since they sum up 
to zero, each of them must be constant, i.e.. 



d 

(4.10) |a,-(x;)fa/'(xO = T,-, £t; = 

(=1 



C'' 3 -regularity of each a, follows from Corollary 13. 2l and the proof of Proposition 14.2 l is concluded. 

□ 

In a number of geometrical and physical problems, it is often that solutions inherit radial sym- 
metry from the data. It is therefore interesting to analyze the regularity theory for radially symmetric 
solutions. 

Proposition 4.3. Let u € C^{B\) be radially symmetric, i.e., u(X) = v{r)for r := and satisfy 

AcxjM = inB\ 

in the viscosity sense. Assume f € C°([0, 1] x M). Then v G C^^c . 
Proof. Direct computations yield 

V«(X) = v'{r)^ 



D\{X) = -^v"{r)X®X + v'{r) 

LA. 



1 1 

—\d-—-X®X 

\x\ 



Therefore, we conclude v satisfies, 

(4.11) v'{rf-v"(r)=f{r,v{r)), 
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in the viscosity sense. Since 

fir,v{r))€LUB,), 

we can apply Corollary l3.2l to assure v € Cj^J^" and the Proposition is proven. □ 
We comment that this regularity from Proposition 14. 3 l is optimal as 

4 

AoojXp = const. 

Let us further mention that for functions with bounded infinity-laplacian, E. Lindgren, following 
ideas from Q, has recently established Lipschitz estimate and everywhere differentiability. 

4.3 Further degenerate elliptic equations 

Another interesting example to visit is the /7-laplacian operator, p>2: 

(4.12) Apu:=div{\Vu\P^^Vu) . 

It appears for instance as the Euler-Lagrangian equation associated to the /^-energy integral 

J {DuYdX min. 

Equations involving the the p-laplacian operator has received a great deal of attention for the past 
fifty years or so. In particular, the regularity theory for p-harmonic functions has been an intense 
subject of investigation, since the late 60's, when Uraltseva in |[26l proved that weak solutions to the 
homogeneous 77-laplacian equation 

(4.13) Aph = 0, 

is locally of class c''"('^''^'\ for some a{d,p) > 0. The sharp regularity for p-harmonic functions 
in the plane was obtained by Iwanec and Manfredi, |[T3l . The precise optimal Holder continuity 
exponent of the gradient of p-harmonic functions in higher dimensions, d >3, has been a major open 
problem since then. 

The /7-laplacian operator can be written in non-divergence form, simply by passing formally the 
derivatives through: 

(4.14) ApU = \V\P-^Au + {p-2)\Vu\P-'^A^u. 

The notions of weak solutions, using its divergence structure in (14. 12K and the non-divergence form 
in (14.141 ) are equivalent, lfT4l . 

Within the context of functions with bounded /^-laplacian, the conjecture is that the optimal reg- 
ularity should be C'^' , where 

P P 

Our next result gives a partial answer to this conjecture. 
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Proposition 4.4. Let p >2 and u satisfy 

|ApM| < C, in B\. 

Assume u is radially symmetric, i.e. u{X) = v{r)for r := Then v G Cy^{B\). 
Proof. Direct computations as in Proposition I4.3l yield 

I l\V—2\ //III /|p— 4 /2| ll\ ^ 

\v Y \^ \ ^\^\ • V |v I < C, 
in the viscosity sense. Thus, an application of Corollary |3.2| reveals that 



and the Proposition is proven. □ 

Estimates of the form |Vm|''^^|D^m| < C ai^e not rare in a number of problems involving the 
/7-laplacian operator, see for instance WT\ . Let us also mention that estimates of the form {e + 
are usually obtained for bounded weak solutions to divergent form equations, 
Di (A'(Dv)) = 0, see for instance yjOi], Chapter 8. 

Let us finish this Section, by revisiting the derivation of the infinity-laplacian operator as the limit 
of p-laplace, as p ^ oo. Let h G C^{B\) be an infinity-harmonic function. For each p ^ 1, let hp be 
the solution to the boundary value problem 

f \php = 0, inB3/4 
y hp = h, on 8^/4. 

It is known that hp form a sequence of equicontinuous functions and hp ^ h locally uniformly to h. 
In particular 

Aco/ip = o(l), asp^oo. 
Hereafter, let us call hp the p-harmonic approximation of the infinity-harmonic function h in B3/4. 

Proposition 4.5. Let h gC^\Bi) be an infinity -harmonic function and hp its p-harmonic approxima- 
tion. Assume \l^^hp \ = 0{p^^) as p ^ 0°. Then h {81/2)- 

Proof. Since hp is p-harmonic, it satisfies 

\VhpfAhp = {2-p)AJip. 

By the maximum principle, 

l|/^nll/~ < • 



From the approximation hypothesis and Corollary 13.21 we deduce 



\hn\\ I 1 < C, 

"'c''MBl/2) - ' 



for a constant C that is independent of p. The proof of Proposition follows by standard reasoning. □ 
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Another interesting Proposition regards /7-harmonic functions with bounded infinity-laplacian. 
Proposition 4.6. Let u be a p-harmonic function in B\ C W^. Assume ^u G L°°{Bi). Then u G 

Proof. The proof follows by similar reasoning as in the proof of Proposition |43] We omit the details. 

□ 



We leave as an open problem whether Proposition 14. 6 1 holds true without the extra assumption 
on the boundedness of the infinity-laplacian. We further conjecture that if a{d,p) is the optimal 
(universal) Holder continuity exponent for p-harmonic functions, then 

o;((i, 7?) > ^ + o(l), asp ->oo. 

5 Universal compactness 

From this Section on we start deUvering the proof for the main sharp regularity estimate announced 
in Section [3l namely, Theorem 13.11 In this first step, we obtain a universal compactness device to 
access the optimal regularity theory for solutions to Equation ( 11.21 ). The proof we shall present here 
uses the main technical tool obtained in the recent work of Imbert and Silvestre, lITTI . 

Lemma 5.1. Let q € be an arbitrary vector and u € C{B\), a viscosity solution to 
(£q) \d, + yu\^F{X,D'u)=f{X), 

satisfying < 1. Given 5 > 0, there exists £ > 0, that depends only upon d,)i,A, and y, such 

that if 

(5.1) ||M|r^|lf(X,M)-F(0,M)||i.(5.) + ||/||L-(B0<^^' 

then we can find a function h, solution to a constant coefficient, homogeneous, {X, A) -uniform elliptic 
equation 

(5.2) 5(£>'/i)=0, Si/2 
such that 

(5-3) ll"-^llL-{ei/2) < 5. 

Proof. Let us suppose, for the sake of contradiction, that the thesis of the Lemma fails. That means 
that we could find a number 5o > and sequences, Fj{X,M), fj, qj and uj, satisfying 

(5.4) Fj{X,M) is (A, A) -elliptic, 

(5.5) ||M|r' • ||F(X,M) -f (0,M)||i.(B,) = o(l), 

(5.6) ll/illL-(B,)=o(l)' 

(5.7) ||"7I|l"(Bi) < 1 and \qj + Vuj\'^Fj{X,D^Uj) =fj, 
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however, 

(5.8) sup \uj — h\ > do, 



'1/2 

for any h satisfying a constant coefficient, homogeneous, (A, A) -uniform elliptic equation (I5.2I ). 

Initially, arguing as in ifTTI . the sequence uj is pre-compact in C°(B 1/2) -topology. In fact, as in 
Ifm . Lemma 4, there is a universally large constant Aq > 0, such that, if for a subsequence {qjj}jteN> 
there holds, 

|qAl>^o, yken, 

then, the corresponding sequence of solutions, {uj^}j(^^, is bounded in C*^'' (^2/3)- If 

|qy|<Ao, yj>jo, 

then, by Harnack inequality, see |[T2l . {My}y>yo is bounded in C*''^ (^2/3) for some universal < j3 < 1. 
From the compactness above mentioned, up to a subsequence, 

Uj — > Moo, locally uniformly in B2/3, 

for some Uoo G €^'^{82/3). Our next goal is to prove that the limiting function Moo is a solution to a 
constant coefficient, homogeneous, (A, A)-uniform elliptic equation (15.21) . For that we also divide 
our analysis in two cases. 

If \qj \ bounded, we can extract a subsequence of {qj}, that converges to some qco S M'^. Also, by 
uniform ellipticity and (I5.5l l. up to a subsequence Fj{X,-) -^^{■), and 

\qoo + Vu^\^^{D^Uc.) = 0, 

Arguing as in ifTTI . Section 6, we conclude that Moo is a solution to a constant coefficient, homogeneous 
elliptic equation, which contradicts (15. St . 

If [q^l is unbounded, then taking a subsequence, if necessary, \qj\ 00. In this case, define 
ej = q;/|qy | and then uj satisfies 



Iqyl 



Letting j ^ °° and taking another subsequence, if necessary, we also end up with a limiting function 
Moo, satisfying 5oo(D^Moo) = for some (A, A)-uniform elliptic operator, Joo. As before, this gives a 
contradiction to (15.81 ). The Lemma is proven. □ 



6 Universal flatness improvement 

In this Section, we deliver the core sharp oscillation decay that will ultimately imply the optimal C' " 
regularity estimate for solutions to Equation (ll.2l i. The first task is a step-one discrete version of the 
aimed optimal regularity estimate. This is the contents of next Lemma. 
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Lemma 6.1. Let q^^W^ be an arbitrary vector and u G C(Si) a normalized, i.e., \u\ < 1, viscosity 
solution to 

(£q) \q + Vu\^FiX,D^u)=f{X). 

Given a G (0, Ob) n (0, ■:^], there exist constants < Po < 1/2 and Eq > 0, depending only upon 
d,X,A,Y and a, such that if 

(6.1) ||M|r^ • \\F{X,M)-F{0,M)\\L-iB,) + ||/||L-{eo < £o, 

then there exists an ajfine function i{X) = a+b ■ X, such that 



sup|m(X)-£(X)| <P(5+" 

Furthermore, 



\a\ + \b\ <C{d,?i,A), 

for a universal constant C{d,X,A) that depends only upon dimension and ellipticity constants. 

Proof. For a 5 > to be chosen a posteriori, let h be a solution to a constant coefficient, homo- 
geneous, (A, A) -uniform elliptic equation that is 5-close to u in L°°(Bi/2)- The existence of such a 
function is the thesis of Lemma [STTl provided Eq is chosen small enough, depending only on d and 
universal parameters. Since our choice for 5 - later in the proof - will depend only upon universal 
parameters, we will conclude that the choice of £o is too universal. 
From normalization of u, it follows that 

||/j||l-(b,/2)<1 + 5<2; 

therefore, from the regularity theory available for h, see for instance [2], Chapters 4 and 5, we can 
estimate 

(6.2) sup|/j(X)-(V/i(0)-X + /i(0))| < C{d,X,A)-?+'^ Vr>0, 

Br 

(6.3) \Vh{0)\ + \h{0)\ < C{d,X,A), 
for a universal constant < C(J, A,A). Let us label 

(6.4) i{X) = Vh{0)-X + h{0). 
It readily follows from triangular inequality that 



(6.5) sup\u{X) - £{X)\ < 5 + C(J,A,A)-Po+'*. 
Now, fixed an exponent a < ao> we select po and 5 as 

(6.6) Po := "0-?/ 

(6.7) 5 : 



2C(J,A,A)' 

l+g 

1 / 1 \ "0-" 



2 \2C{d,X,A) 
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where < C(<i, A, A) is the universal constant appearing in (I6.2I ). We highlight that the above choices 
depend only upon d,X,A and the fixed exponent < a < tto- Finally, combining (I6.2I ). (I6.5I ). (16.61 ) 
and (16.71) . we obtain 



2 ^ s l+« 



sup|«(X)-£(X)| < -(^^^^j«-''+C(J,A,A).po^+«.po" 



Po 

and the Lemma is proven. □ 

In the sequel, we shall iterate Lemma [6ni in appropriate dyadic balls as to obtain the precise sharp 
oscillation decay of the difference between u and affine functions i^. 

Lemma 6.2. Under the conditions of the previous lemma, there exists a sequence of affine functions 
ik{X) : = ak + bk-X satisfying 

(6.8) \ak+i-ak\+p^o\h+i - < CoP^'+"^^ 

such that 



(6.9) sup|KX)-4(^)|<Po^'+"^ 

B k 
Po 

where a is a fixed exponent within the range 

(6.10) aG(0,ao)n(o, ^ 



1 + 7. 

and Co is a universal constant that depends only on dimension and ellipticity. 

Proof. We argue by finite induction. The case k=lis precisely the statement of Lemma |6.1| Suppose 
we have verified ( 16.91 ) for 7 = 1,2, •• • ,k. Define the rescaled function 

Po 

It readily follows from the induction assumption that |v| < 1. Furthermore, v satisfies 

\p-'"'h + yv\^ Ft{X,D\) = f,{X), 



where 



and 



Mx) = P^^'-''''^'^^fiP^x) 



F,{X,M) := pJ('-«V i^p^X, -j^M 
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It is standai^d to verify that the operator f/t is (A, A)-elliptic. Also, the OD-norm of the corresponding 
coefficient oscillation of Fk, as defined on (I2.4I ). hereafter called jS/t, does not increase. Also, one 
easily estimate 

(6.11) < pf-''^'^'Hf\\L^iB 

Po 

Due to the sharpness of the exponent selection made in (16.101 ). namely a < j^, we conclude {Ft,fk) 
satisfies the smallness assumption (16. IK from Lemma l6n 

We have shown that v is under the hypotheses of Lemma ISTl which assures the existence in the 
affine function i(X) ■.= a + b-X with \a\ + \b\ < C{d,X,A), such that 

(6.12) sup|v(X) < p^+". 

In the sequel, we define the (^+ l)th approximating affine function, {X) := at+i +bk+i X, where 
the coefficients are given by 

ak+i:=ak + plj^^"^''a and bu+i ■=bk + Pq^I. 

Rescaling estimate (16.121 ) back, we obtain 



sup iM(x)-4+i(x)i<pr"^'^"^ 

and the proof of Lemma l6^ is complete. □ 

7 Smallness regime 

In this Section we comment on the scaling features of the equation that allow us to reduce the proof 
of Theorem 13. II to the hypotheses of Lemma |6TT] and Lemma [6^ 
Let V G C(Bi) be a viscosity solution to 

jr{X,Vv)F{X,D^v)=f{X), 

where Jif satisfies (11.31 ) and f is a (A, A) -elliptic operator with continuous coefficients, i.e., satisfying 
(12.11 ) and (12.31 ). Fix a point Yq E and define w : Si ^ M as 

for parameters rj and T to be determined. We readily check that u solves (in the viscosity sense) 

J^^4X,Vu)F^,,{X,D\) = fr^^riX), 
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where 

(7.1) fr,,T(X,M) 

(7.2) =^r,,r(^,p) 
(V.3) fr^,,{X) 

Easily one verifies that F,j ^ is uniformly elliptic with the same ellipticity constants as the original 
operator F, i.e, it is another (A, A) -elliptic operator. Also J^,t satisfies the degeneracy condition 
(11.31) . with the same constants. Let us initially choose 

T:=max{l,||v||z,=o(B,)}, 

thus, [m| < 1 in Bi. Now, for the universal £o appearing in the statement of Lemma l6.11 we select 

r7:=min{l,A-(£o||/||,i)^,«-' (|)}. 

With these choices, u is under the assumptions of Lemma l6. II Proving universal C^-^ regularity for u 
at gives the appropriate C^''^ estimate for the original function v in terms of ||v||/,»(Bj) and ||/||z,-(b,) 
at the generic point Fq £ 5i/2- 

The above reasoning certifies that in order to show Theorem 13. 11 it is enough to work under the 
smallness regime requested in the statement of Lemma 16.11 Once established the desired optimal 
regularity estimate the normalized function u, the corresponding estimate for v follows readily. 

8 Sharp local regularity 

In this Section we conclude the proof of Theorem 13. II From the conclusions delivered in Section|7] it 
suffices to show the aimed C^-^ estimate at the origin for a solution u under the hypotheses of Lemma 
16. H and Lemma |6^ For a fixed exponent a satisfying the sharp condition (16.101 ). we will establish 
the existence of an affine function 

such that 

and 

sup|m(X) -4(^)1 <Cr'+«, Vr<l, 

Br 

for a constant C that depends only on J, A, A, 7 and a. 



r,r+2 
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Initially, we notice that it follows from (16.81 ) that the coefficients of the sequence of affine func- 
tions £k generated in Lemma \62\ namely bk and at, are Cauchy sequences in and in M, respec- 
tively. Let b^, and be the limiting coefficients, i.e.. 



(8.1) 
(8.2) 



lim bic 
lim Qk 



It also follows from the estimate obtained in (16.81 ) that 

(8.3) < 

(8.4) \k-bk\ < 



_Q>_k{l+a) 



Now, fixed < r < po, we choose ^ G N such that 



We finally estimate 



sup|m(X)-4(X)| < sup 4(^)1 

Br 



B k 



< sup \u{X) - 4(X) I + sup \hiX) - L{X)\ 



B k 
Po 



B k 
Po 



< 



< 



k{l+a) , Co m+a) 

Po — -Po 



1 



Po 



1+a 



1 + 



1-Po 
Co 



1-Po 



J+a 



and the proof of Theorem l3.1l is complete. 



□ 
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